Non-invasive characterization of human vasculature

ABSTRACT

Vascular conditions are detected non-invasively in the human body using a collection of information from small local regions of the vasculature. An array of accelerometers or other sensors are attached to the head or other points of interest of a patient and blood flow sounds are recorded. Vibration signatures of vessel structures such as branches, aneurysms, stenosis, etc. using random, periodic, band limited or transient analysis provides a library for further processing. The signature library is used to localize the origin of the recognized vascular feature, and the localized feature is presented to the physician in a clinically relevant manner.

This application claims benefit of provisional application Ser. No. 60/838,624, filed Aug. 17, 2006.

BACKGROUND OF THE INVENTION

This invention concerns detection of conditions of human vasculature non-invasively, from outside the body, and in particular, locating aneurysms, partial or complete stenoses, ruptures, peripheral bleeding and other abnormal conditions pertaining to the cerebral vasculature, profiling blood or other fluid flow through or around the cerebral vasculature and mapping the resultant flow to provide a spatially resolved characterization of the major vasculature within the head.

The system and method of the invention are also useful to detect vascular conditions in other parts of the body, especially the limbs and areas of other major vessels. Stroke is a manifestation of vascular injury to the brain which is commonly secondary to atherosclerosis or hypertension, and is the third leading cause of death in the United States. Stroke can be categorized into two types, ischemic stroke and hemorrhagic stroke. Additionally, a patient may experience transient ischemic attacks, which are in turn a high risk factor for the future development of a more severe episode.

Examples of ischemic stroke encompass thrombotic, embolic, lacunar and hypoperfusion types of strokes. Thrombi are occlusions of the arteries created in situ within the brain, while emboli are occlusions caused by material from a distant source, such as the heart and major vessels, often dislodged due to myocardial infarct or atrial fibrillation or carotid disease or surgical or percutaneous intervention. Thrombi or emboli can result from atherosclerosis or other disorders, for example, arteritis and lead to physical obstruction of arterial blood supply to the brain. Lacunar stroke refers to an infarct within non-cortical regions of the brain. Hypoperfusion embodies diffuse injury caused by non-localized cerebral ischemia secondary to low cerebral perfusion, typically caused by myocardial infarction, arrhythmia, blood loss or prolonged low blood pressure.

Hemorrhagic stroke is caused by intra cerebral or subarachnoid hemorrhage, i.e., bleeding in to brain tissue, following blood vessel rapture within the brain or venous thrombosis. Intra cerebral and subarachnoid hemorrhages are subsets of a broader category of hemorrhage referred to as intracranial hemorrhage. Intra cerebral hemorrhage is typical due to chronic hypertension and a resulting rupture of an arteriosclerotic vessel. Stroke-associated symptom(s) of intra cerebral hemorrhage are abrupt, with the onset of headache and steadily increasing neurological deficits. Nausea, vomiting, delirium, paralysis, seizures and loss of consciousness are additional common stroke-associated symptoms.

In contrast, most subarachnoid hemorrhage is caused by head trauma or aneurysm rupture which is accompanied by high pressure blood release which also causes direct cellular trauma. Prior to rupture, aneurysms may be asymptomatic, or occasionally associated with headaches. However, headache typical becomes acute and severe upon rupture and may be accompanied by varying degrees of neurological deficit, vomiting, dizziness, and altered pulse and respiratory rates, phobophobia and severe headache and or neck stiffness.

Current diagnostic methods for stroke include costly and time-consuming procedures such as non-contrast computed tomography (CT) scans, electrocardiogram, magnetic resonance imaging (MRI) and angiography. Determining the immediate cause of stroke is difficult. CT scans can detect parenchymal bleeding greater than 5 mm and 95% of all subarachnoid hemorrhages. CT scans often cannot detect ischemic strokes until 6 hours from onset, depending on infarct size. CT only identifies 48% of acute strokes in the first 48 hours. MRI may be more effective than CT scan in early detection of ischemic from hemorrhagic stroke, and is not widely available. Angiography is a definitive test to identify stenosis or occlusion of large and small cranial blood vessels, and can locate the cause of subarachnoid hemorrhages, define aneurysms, and detect cerebral vasospasm. It is, however, an invasive procedure and is also limited by cost and availability.

Immediate diagnosis and care of patient experiencing stroke can be critical. For example, tissue plasminogen activator (tPA) given within three hours of symptom onset in ischemic stroke is beneficial for selected acute stroke patients. In contrast, thrombolytics and anticoagulants are strongly contraindicated in hemorrhagic strokes. Thus early differentiation of ischemic events from hemorrhagic events is imperative. Moreover, delays in confirmation of stoke diagnosis and identification of stroke type limit the number of patients that may benefit from early intervention therapy. In addition, continuous monitoring of stroke patients is not possible with CT or MRI scanners thus only one snapshot is time is available for diagnosis and treatment. Clinical observations are the basic tool that is used to monitor the progress of stroke patients.

Early detection of an aneurysm is beneficial as it can frequently be treated either by surgical procedure of clip occlusion or by endovascular coil embolism. Presently, approximately three quarters of patients are treated with clip occlusion, the remainder with endosvascular coil embolism. Either surgery, particularly the endovascular procedure, can by performed with low complication rate and high rate of success.

Once an aneurysm ruptures, however, the patient declines rapidly due to major brain injury, and over 50% of aneurysm rupture patients die acutely. Thus detection of at risk aneurysms is of great benefit. Thus, a physician faced with possible aneurysm warning signs must judge whether the symptoms warrant the trauma, expense and morbidity of contrast angiography.

Ferguson, in J. Neurosurg. 36:560-563 (1972), suggested detecting aneurismal signals by recording sounds from aneurysms exposed at operations using a cardiac phono-catheter.

Kosugi et al., in Stroke 14 (1) 37-42 (1983), disclosed the use of a “cement wall microphone”, (contact accelerometer) in contact with the cranium and the teeth in an attempt to detect aneurysms.

Despite these known devices, there remains a need for an acoustic detector that is designed to more effectively record, analyze, localize and present, in a clinically useful manner, cerebral arterial and venous conditions. Devices and methods to date have not localized the detected vascular conditions and do not present the collected data in a clinically useful way.

It is the purpose of this invention to provide a noninvasive method for differentiating ischemic from hemorrhagic stroke and to continuously monitor the condition of stroke patients by providing clinically useful localized information. It is also the purpose to monitor or screen high risk patients to detect conditions that may lead to a first or subsequent stroke.

SUMMARY OF THE INVENTION

The current invention is a system and method for detecting a vascular condition non-invasively in the human body, by a collection of signal information from local, small regions of the vasculature. This is accomplished in a preferred embodiment by attaching, or contacting an array of accelerometers, or other sensors, to the head, and/or other points of interest, to a patient and recording signals. The vibration signatures of blood vessel structures such as branches, aneurysms, stenosis or other structures using random, periodic, band limited or transient analysis provides a library for further processing. The library becomes part of a cerebral modeling arithmetic computer using basis functions, or other artificial neural network. The developed signature library is then used to facilitate localizing the origin of the recognized vascular feature; the localized feature is then presented to the physician in a clinically relevant manner.

The data to build the neural network (NN) and/or data store for correlation is obtained by acquiring data from patients with known pathologies as well as controls, and optionally with studies of constructed or cadaver vasculatures. The data is formatted according to the requirements of each processing method before being used for training or algorithm construction.

These and other objects, advantages and features of the invention will be apparent from the following description of a preferred embodiment, considered along with the accompanying drawings.

DESCRIPTION OF THE DRAWINGS

FIG. 1A is a plan and elevation view of sensors on a head.

FIG. 1B is a block diagram of the steps of the invention.

FIG. 2 is a block diagram of the sensor processing steps.

FIG. 3 is a diagram of the correlation circuit.

FIG. 4 is a diagram of a trigger pulse amplitude VS time

FIG. 5 is a diagram of sensors at various distances from a pulse source.

FIG. 6 is a diagram of a trigger pulse and a received pulse at a distant sensor.

FIG. 7 is the overlay of the trigger pulse and the received pulse shifted by time delta t1.

FIG. 8 is a diagram of a trigger pulse and a received pulse from a further sensor.

FIG. 9 is a diagram of a trigger pulse and a received pulse from the furthest sensor.

FIG. 10 is a table the time of first signal arrival from the trigger pulse for each sensor.

FIG. 11 is a diagram of two intersecting spheres of different diameters.

FIG. 12 is a diagram of three sensor signals shifted according to the table in FIG. 10.

FIG. 13 is a diagram of the sum of the three sensor signals in FIG. 12.

FIG. 14 is a received sensor signal as shown in FIG. 6 and typical noise signal.

FIG. 15 is a diagram of the sum of the noise signal and the signal from sensor 1.

FIG. 16 is a diagram of the sum of noise signal and the signal from sensor 2.

FIG. 17 is a diagram of the superposition of the noise signals from each of the three sensors.

FIG. 18 is a diagram of the sum of the noise signals from each of the three sensors.

FIG. 19A is a plan view in one plane of four sensors on a human head with the distances from a noise source to each sensor.

FIG. 19B is a plan view in one plane of four sensors on a human head with several possible sources of sound indicated.

FIG. 20 is a plan view of the head with sensors shown and the geometry of beam forming.

FIG. 21 is a table giving the distances for four sensors and five beam forming position.

FIG. 22 is a plot of the width of a blood vessel determined by beam forming.

FIG. 23 is a diagram indicating the width of a blood vessel based on merging the NN output and the beam forming data.

FIG. 24 is a figure of a patient being injected with ultrasound contrast agent.

FIGS. 25A and 25B are figures of a branch in a blood vessel without and with ultrasound contrast.

FIG. 26 is a figure of an aneurysm in the branch of a vessel as imaged with bursting bubbles.

FIG. 27 is a figure of what an image of the vessel and aneurysm would be displayed after processing.

FIG. 28 is a figure of a branching vessel and areas of perfusion and non-perfusion.

FIG. 29 is a figure of a vessel blocked by a stenosis.

FIG. 30 is a figure of a AVM (arterial venous malformation).

FIG. 31 is a diagram of how bursting bubbles combine with a NN image to produce an improved image.

FIG. 32 is a figure that depicts how two segments derived from imaging can be connected by knowledge of the anatomy.

FIG. 33 is a figure that depicts how two segments derived from imaging can be connected by bursting bubble data.

FIG. 34 is a figure that depicts an area of pooled blood in proximity to an aneurysm.

FIG. 35 is a flow diagram outlining the overall system.

DESCRIPTION OF PREFERRED EMBODIMENTS

General Outline of Process

Periodic pulses from the heart produce waves of expanding blood vessels within the body. The tissue surrounding the blood vessels is displaced during the pulse and contracts again after the pulse. This displacement propagates outward from the blood vessel. In the case of the brain, this displacement reaches the skull and displaces the bone in response the displacement of the blood vessel wall. Extremely sensitive accelerometers record this displacement. The character of the signal recorded by the accelerometers is dependent on the nature of the displacement caused by the blood pulse. At a restriction in the vessel the displacement before the restriction is larger than it would be without the restriction and the displacement beyond the restriction is less. The spatial distribution of displacement produces a different signature than that recorded by an unrestricted vessel. Likewise an aneurysm allows a circulation of blood within the bulb of the aneurysm during the pulse and produces a periodic signal that modifies the displacement signal. Other geometric arrangements of vessels likewise produce unique signatures in displacement.

The accelerometers are very sensitive, typically 500 mV/g or more. Such an accelerometer is a model from Dytran, Calif.

Limitations of Acoustic Signal Detectors

Accelerometers have a distinct advantage over microphones since microphones require acoustic signals to propagate through the soft tissue of the brain, conduct through the hard bone of the skull, again transfer to the soft tissue of the skin and finally be transferred to the microphone. Such devices have been most successful when attached to the teeth, positioned within the ear or focused on the eye sockets. Thus the use of microphones has been more successful when they can either have direct bone contact or look where there is no skull present. Accelerometers measure the displacement of the entire structure.

Elimination of External Sources of Signal

The microphone is also very sensitive to sound signals produced from outside the body, such as sirens, talking, doors slamming and other sound signals. Accelerometers are much less sensitive to acoustic signals and generally are only sensitive to signals along the primary axis of the detector.

Localization with Phase

The use of multiple accelerometers positioned around the head allows signals from different parts of the brain to be distinguished by location. Beyond just the relative magnitude of the signals from each of the sensors the phase relationship of the signals between sensors is used to determine the location of the displacement source. Since there are a multitude of signals produced during the pulse from different spatial positions and with the pulse traveling from the proximal to distal portions of the vasculature, a means to distinguish each of these signals is needed. This differentiation is accomplished by using a priori information about the head and the approximate positions of each of the accelerometers. For example signals can not emanate from locations outside of the head nor at times that are not related to the blood pulse.

Timing for Localization

A gating sensor, an EKG electrode, is used to provide timing information to the system. Since the heart rate varies in a given person both at rest and with exertion, the timing signal is used as an approximate starting point for timing the signal arrival at the sensors but cannot determine the pulse timing. Leeway is still needed as the intra pulse timing may be increasing or decreasing for each pulse.

Localization Example

Signals data is captured from multiple sensors during a pulse. The signals from the multiple sensors are shifted in time by a guessed amount and then added together. The resulting signal is then analyzed either by direct amplitude, by Fourier transform or by another algorithms. This process is continued over the valid range of phases from each sensor. If a given signal is correlated with the signal from another sensor then the resulting combined signal will have an enhanced signature when the phases represent the location of the signal and not for surrounding phases within the allowable phase space. The quality of a signature can be quantified and that quantification can be established for each phase shift between sensors. As the phase is shifted over the possible range the figure of merit (quantified result of a single phase relationship) will increase and decrease. The peaks represent recognized signatures at localized positions. While this process is straightforward it is very processor intensive. The allowable phase shifts for all the sensors can be calculated beforehand and this range used to limit the correlation space. The computational problem can be parsed such that parallel processors tackle the problem in parallel thus reducing the time from data collection to data presentation. It is important to note that once the initial data has been analyzed the relative phase between sensor data for the many signals within the brain is known, continued analysis of captured data is much less computational intensive. This a priori knowledge vastly reduces the computational requirements of the system and allows near real time presentation of the data. Current computers are certainly capable of completing the vast amount of computations in times suitable for clinical use and this will just improve as the speed and processing throughput of computers increases.

Creation of Synthetic Image

The data needs to be presented in a recognized manner for the clinician to make use of it. By the nature of the data capture and processing, the most prominent structures are most readily identified and localized. This is not a complete map of the brain vasculature. Using a priori knowledge of the typical structure of the vasculature (and noting that there are significant deviations from the “typical” structure) the identified structures can be placed onto a cartoon of the entire vasculature, replacing the typical cartoon representations with improved representations based on the data analysis. This synthetic image is then presented to the clinician; as a 3D representation, a simulated CT or MRI scan or as a angiogram. As data continues to be acquired and processed the synthetic image can be updated.

Ultrasound Bubbles

Ultrasound contrast agents are materials that can be injected as an IV to improve the imaging of blood vessels by ultrasound. The contrast consists of very fine particles and a method of attaching very small bubbles to the particle. To avoid any associated problems with embolisms the particle and bubbles together are smaller than a typical red blood cell thus allowing the bubbles to pass through from the arteries to the veins. These bubbles are also engineered to not join together as they are formed about a tiny particle rather than free standing bubbles that might have a tendency to enlarge by combining. The purpose of injecting bubbles into the vasculature is to increase the ultrasonic contrast of the blood by adding well dispersed low density material. Contrast agent is rapidly eliminated in the blood by the breaking of the newly formed bubbles.

An unrecognized feature of the contrast agents is that as each bubble breaks it creates a disturbance within the vasculature which is highly localized. Ultrasound contrast media are already FDA cleared for use in stroke diagnostics and vessel imaging. To ensure the safety of injecting bubbles into the blood the bubble size is very well controlled. Bubbles of the same size have a well defined and distinctive signature when they break. This consistency allows for precise localization of the resulting signal. Since the bubbles are smaller than a red blood cell, the bubbles are able to perfuse from the arteries to the veins. The bubbles break randomly over a controlled period of time, typically five to 20 minutes after injection. If the amount of contrast media that is injected is reduced significantly from the amount typically given for ultrasound imaging the number of bubbles that break can be controlled to have, on average, a few milliseconds between bursts. Thus a system that is sensing the breaks from multiple sensors can localize the position of the bubble when it broke to a few mm or better resolution. Summing over a large number of events will produce a 3D map of where the bubbles have reached. One of the major objectives of stroke treatment is determining whether a patient has a stroke at all and if they have had a stroke to differentiate an ischemic from a hemorrhagic. This differentiation is critical to caring for the patient. If clot dissolving drugs are used on a hemorrhagic patient the likely outcome is significantly worse, including death, over no treatment at all. Likewise not treating an ischemic stroke patient within the first three hours (for current drugs) of the stroke results in no treatment at all.

A 3D bubble map will show areas that are not perfused, an area with little or no bubbles showing up; the signature of a vessel blockage or ischemic stroke. Aneurysms, the typical cause of hemorrhagic strokes, will show up in both the direct sensor signals and in the 3D bubble map. A lack of perfusion in an area around the aneurysms would be a strong indicator of non flowing pooled blood. A perfused area distal to the aneurysm would be a further indication of no ischemic stroke. Thus the system would allow rapid, clinically relevant, differentiation of stroke type enabling treatment to be delivered in a timely manner.

Merger of Data

The sensors deliver data on major unusual events in the vasculature. Merging this data into a “typical” 3D map of the complete vasculature ignores the real differences between individuals. It still provides clinically relevant data, such as lack of blood flow in an area of the brain and then a subsequently reperfusion of that area. This highlights the major advantage of having a continuous monitoring system over a snapshot monitoring such as CT, MRI or contrast angiography. Adding the bubble map would allow a complete map of the vasculature since the bubbles burst through out all the perfused vessels within the brain. Since the same sensors, and sensor locations, are used to produce both maps, the two data sets can be merged to form a single complete map of the vasculature. In addition, this data set can be merged with the CT scan or other scans to combine the data from each to provide the clinician a better picture of the patients condition.

Neural Network

The array of different accelerometer signatures produced by features within the vasculature is quite large due to the variable physiology of patients but there are generally only a few major underlying features. This type of data is ideal for using neural networks to identify the features causing the signature to be categorized. As the library of unique features and the associated signatures grows the neural network will improve in correctly identifying features in patients. Neural networks typically are trained by inputting a set of known inputs and known outputs and allowing the weights of the neural connections to change to optimize the matching of the inputs and outputs. When a new input is presented to the neural network the output closest to the input is given the highest output even if the input is not a perfect match to any of the training set.

Beam Forming

In a process similar to localization previously discussed, beam forming varies the phase between data produced by different sensors. It is the purpose of beam forming to systematically vary the phase and retain the resulting signal. With small steps in the phase a very high spatial resolution map of the vasculature is produced. Beam forming is used when a feature of interest is localized and more detailed information is desired about that feature. Beam forming is computational intensive but lends itself to parallel processing and will be aided by improvements in processing power of both general purpose and special purpose processors.

In the drawings, FIG. 1A shows a patient with sensors attached to various places on the head. The sensors are attached by cables, but could be connected wirelessly, to the analysis portion of the system. FIG. 1B shows a block diagram of the analysis and visualization portion of the system. The signal arrives from each sensor and is conditioned and processed in the data processing block. The processed signal is digitized and processed in parallel through each of the four signal type processors using the causal neural net. These processors analyze the signal for matches to library conditions and identify matches. The matched signals from multiple sensors are then processed by the signal localizer to localize the source of the matched signal. This process is repeated for each identified signal both by signal type and for multiple locations. Patient data is input into the control processor and provided to the causal neural network to further improve identification and localization. The collection of all the localized and identified data is processed by the image database for presentation. The digitized input signal is also stored in the array data memory for further analysis. Localized features can be further characterized by using beam forming techniques. In the beamforming mode multiple sensors are used to “look” at the same location and the resulting signals analyzed for additional information. Beam forming can be used to shift or scan the region of interest to further characterized and localize the feature.

Some of the features of the system as represented in FIG. 35 are:

-   -   The use of noise cancellation.     -   The elimination of overriding signals.     -   The localizers for each signal type.     -   The use of the raw data memory.     -   The causal network.     -   The patient factors input and how it is used.     -   The image database.     -   The image.     -   The dynamic filtering and its purpose.

FIG. 2 is a more detailed block of the sensor and data acquisition block in FIGS. 1A and 1B. The sensor signal is amplified and possibly filtered, sometimes dynamically, before being amplified and finally digitized. The signal from each of the sensors in FIG. 1A are digitized maintaining the time relationship between each sensor.

FIG. 3 shows a block diagram of the correlator bank that localizes the signals. All the sensor signals, with their timing information, are fed into to the correlator bank. The correlator determines, by changing the time relationship between different sensor signals, when two signals are from the same source and record those time relationships.

Multipath rejection via modeling of the signal phase delay of the skull (or body) may be used to improve accuracy of the algorithms.

The relative attenuation of the intervening materials may also be used.

The signal that is common to the multiple sensor signals is recorded as the signal that was created at the location that is determined by the time differences between the arrival of the signal at each sensor, as will be further explained below. FIG. 4 shows a simple trigger pulse test source signal in amplitude versus time.

FIG. 5 shows a diagram of the location of a signal source and three sensors, S1, S2 and S3, shown at different distances from the signal source.

FIG. 6 shows a time VS amplitude waveform of the trigger pulse and the arrival of the signal at sensor S1. The time that the signal takes to arrive at sensor S1 is delta t1.

FIG. 7 shows the superposition of the two signals in FIG. 6 translated in time by delta t1.

FIG. 8 shows the same data as FIG. 6 but for sensor S2 that is further from the source than sensor S1. The time for the signal to arrive from the source at sensor S2 is delta t2.

FIG. 9 shows that same data as FIGS. 6 and 8 but for sensor S3 with a transit time of delta t3.

FIG. 10 shows a table with the arrival time for the signal for each of the sensors, S1, S2 and S3.

FIG. 11 shows two different radius intersecting spheres with the dotted line showing the locus of points where the distance between the centers of sphere 1 and sphere 2 are constant. This locus of points represents the possible position of the source based on only the signal from two sensors. Adding a third sensor, not shown, would restrict the location of the source to a single location, thus fixing its location with respect to the three sensors. FIG. 12 shows all the sensor signals from FIGS. 6, 8 and 9 translated by the times listed in FIG. 10.

FIG. 13 shows the sum of the three signals shown in FIG. 12.

FIG. 14 shows the typical signal noise that might be expected to be present on the sensor signal along with the sensor signal as shown in FIG. 6.

FIG. 15 shows how the sensor signal in FIG. 14 would look for a typical signal noise environment.

FIG. 16 shows the same signal as FIG. 15 but for sensor S2 rather than sensor S1.

FIG. 17 shows that superposition of the noisy signals from all the sensor, translated by the times shown in the table in FIG. 10.

FIG. 18 shows the sum of all the signals in FIG. 17 with the noise reduced due to the summing of random noise and the signal increased by the superposition of the common signal.

FIG. 19A shows a pictorial diagram of a patient with sensors S₁ through S₄ attached to the head. A signal source at X is located at distances D₁ through D₄ from the respective sensors. A waveform is shown at X depicting the acoustic signature that is emitted at location X.

FIG. 19B shows a pictorial diagram of a patient as in FIG. 19A. Four locations within the head of the patient are shown; points A, B C and C′. A centerline of symmetry is shown with point C′ laying on the centerline and equal distance from sensors S₂ and S₄. Point C is also on the centerline and a short distance from point C′. Small changes in the distances D₂ and D₃ represent the small distance from C′ to C. Beamforming techniques depend on changing the phase relationship between data collected by each sensor in a precise and controlled way, such as represented by the example in FIG. 19B.

FIG. 20 shows a pictorial diagram of a patient with sensors S₁ through S₄ as is shown in FIG. 19A with the additional points A-E shown. Points A through E represent small steps across a point of interest, as depicted by a segment of vessel shown across the points A-E. The distances from each point A-E to each sensor S1-S4 is different and is shown in the table in FIG. 21.

An example would be when giving tPA to dissolve a clot. The clot region could be monitored and tPA administered while watching the blood flow and perfusion so that the right amount of tPA is administered. This would be possible even if the patient was not able to indicate physical improvements, i.e., asleep or in a natural or induced coma. The preferred embodiments described above are intended to illustrate the principles of the invention, but not to limit its scope. Other embodiments and variations to these preferred embodiments will be apparent to those skilled in the art and may be made without departing from the spirit and scope of the invention as defined in the following claims. 

We claim:
 1. A method for detecting, recording, analyzing, and storing data of cranial vasculature of a patient with or without a cranial vascular abnormality, at a region of interest, including the steps of: (a) using an array of accelerometers engaged externally at a plurality of different positions against a patient's head, measuring acceleration and recording pressure wave/motion stimulation acceleration signal data from displacements of the skull at said plurality of positions on the skull, such displacements occurring from heartbeat-induced repeated pulsing of blood flow into the vasculature of the brain and resultant expansion and subsequent contraction of the blood vessels which moves the brain tissue due to such pulsing, the brain tissue movement causing resultant displacement of the skull at said plurality of positions on the skull, and (b) associating said signal data to characteristics of known signals of cranial vascular anomalies, in a computer, and finding any matches of said signal data to known signals of vascular anomalies, and further including (c) correlating signal data from multiple accelerometers of said array to determine a waveform signal common to the signal data from the multiple accelerometers and determining phase time differences among the signal data from the multiple accelerometers to thereby determine origins of a series of signal data collected by the multiple accelerometers from said displacements of the skull, thus localizing at least one signal from the signal data of the multiple accelerometers to define location of a structure of a vascular anomaly, and (d) storing results of the correlation and localizing steps as metadata in a database.
 2. A method for detecting, recording, analyzing, and storing data of cranial vasculature of a patient with or without a cranial vascular abnormality, at a region of interest, including the steps of: (a) using an array of accelerometers engaged externally at a plurality of different positions against a patient's head, the accelerometers having a sensitivity of about 500 mV/g or greater, measuring acceleration and recording pressure wave/motion stimulation acceleration signal data from displacements of the skull at said plurality of positions on the skull, such displacements occurring from heartbeat-induced repeated pulsing of blood flow into the vasculature of the brain and resultant expansion and subsequent contraction of the blood vessels which moves the brain tissue due to such pulsing, the brain tissue movement causing resultant displacement of the skull at said plurality of positions on the skull, and (b) associating said signal data to characteristics of known signals of cranial vascular anomalies, in a computer, and finding any matches of said signal data to known signals of vascular anomalies. 